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ION  BEAM  NEUTRALIZATION  PROCESSES  FOR  ELECTRIC  MICROPROPULSION 

APPLICATIONS 

Adrian  Wheelock*  and  David  L.  Cooke1" 

Air  Force  Research  Laboratory ,  Space  Vehicles  Directorate,  Hanscom  AFB ,  MA 

Nikolaos  A.  Gatsonis* 

Mechanical  Engineering  Department,  Worcester  Polytechnic  Institute,  Worcester,  MA  01609 

Abstract 

While  it  is  common  knowledge  that  ion  beams  are  easily  neutralized  using  a  variety  of  means,  the  precise  process  of 
neutralization  remains  unknown.  With  the  increasing  popularity  of  electric  propulsion,  and  in  particular 
micropropulsion  systems,  this  question  is  of  significant  importance.  Additionally,  it  has  bearing  on  thruster  design, 
space  instrument  calibration,  electrodynamic  tethers,  and  ionospheric  research.  A  review  of  the  present  state  of 
knowledge  on  this  topic  is  presented  as  well  as  results  from  ion  beam  simulations  using  2D  and  3D  ParticIe-in-Cell 
codes.  The  grid  generation  methodology,  adaptation,  charged-particle  transport,  and  field  solver  methodologies  of 
the  3D  code  are  reviewed.  The  simulations  show  electrons  moving  to  neutralize  the  ion  beam  from  background  and 
neutralizer  sources.  The  simulations  show  the  dependence  of  neutralization  on  beam  energy  and  the  electron/ion 
velocity  ratio.  The  results  are  compared  favorably  with  previous  computations  and  experimental  observations. 

flux,  Fi  =  Nft ,  and  the  net  electron  flux, 

Ft  -  X  Nvtih ,  where  N  is  density,  v  is  velocity,  i  and 

e  are  ion  and  electron  subscripts  and  eth  designates 
the  electron  thermal  velocity  for  an  idealized  electron 
source.  Equal  density  and  flux  requires  ve(h  =  4v. .  A 
1  keV  Xenon  beam  has  ^=38,000  m/s  so  a 
matching  electron  velocity  requires  a  source 
temperature  of  about  0.05  eV.  A  challenging,  but  not 
impossible  number,  but  a  collisionless  analysis 
suggests  that  detailed  balancing  is  required,  whereas 
real  systems  quite  easily  achieve  ‘beam  coupling.’ 
Of  course  a  higher  temperature,  lower  density 
electron  source  will  lead  to  a  positive  potential  well 
that  does  trap  electrons,  but  then  the  theory  must 
explain  by  what  process  the  trapped  electrons  shed 
energy  so  as  to  actually  fill  the  well.  Another 
observation  is  that  when  ion  beams  and  neutralizers 
are  operated  in  conducting  vacuum  tanks,  the 
currents  are  closely  coupled  even  though  the 
grounding  tank  eliminates  the  charge  accumulation 
that  could  provide  feedback  for  current  balance  so  it 
appears  that  one  or  more  plasma  mechanisms  must  be 
responsible  for  this  collective  phenomena  ~  charge 
and  current  neutrality  -  which  we  hereafter  call  beam 
coupling. 

Although  the  propensity  for  beam  coupling  is 
clearly  serendipity  for  the  application  of  electric 
thrusters,  our  review  of  the  literature  indicates  it  is 
still  a  conundrum  for  plasma  science.  Further,  new 
electric  micropropulsion  devices  such  as  the  FEEP  or 
the  colloidal  thrusters  or  large  arrays  of  ion  and  Hall 
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INTRODUCTION 

It  is  well  known  that  electrons  injected  into  an  ion 
beam  can  readily  neutralize  it.  The  method  by  which 
the  electrons  are  entrained  and  remain  in  the  ion 
beam  to  neutralize  the  space  charge,  however, 
remains  uncertain.  Proper  modeling  of  the  current 
coupling  and  neutralization  phenomena  will  enable 
development  of  low-current  neutralizers, 
optimization  of  neutralizers  for  electric  propulsion 
devices  as  well  as  aid  in  the  development  of 
emerging  micropropulsion  devices. 

In  the  early  years  of  electric  propulsion  research, 
the  neutralization  question  was  one  of  the 
fundamental  issues  for  successful  development  of 
this  promising  technology.  A  dense  ion  beam 
requires  space  charge  neutralization  to  avoid  a 
potential  barrier  that  can  divert  or  reflect  the  beam. 
The  vehicle  on  which  the  thruster  operates  needs 
current  neutrality  to  avoid  excessive  charging.  In  the 
context  of  collisionless  plasma  theory,  achieving  both 
current  and  charge  neutrality  with  the  same  source  of 
electrons  appears  to  be  nearly  impossible  owing 
mostly  to  the  large  difference  in  mass  between 
electrons  and  the  ions.  For  example,  define  the  ion 
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thrusters  are  still  not  guaranteed  to  behave.  We  might 
also  desire  a  means  to  predict  and  optimize 
neutralizer  operations.  Thus,  a  simulation  technique 
exhibiting  beam  coupling  is  needed.  Additionally, 
results  from  ion  beam  neutralization  modeling  will  be 
applicable  to  ion  beams  for  instrument  calibration, 
electrodynamic  tethers,  ionospheric  research,  and 
fundamental  plasma  physics. 

Our  immediate  goal  is  to  determine  if  what 
might  be  considered  standard  Particle-In-Cell,  PIC, 
techniques  are  adequate  or  if  additional  treatment  is 
needed  to  understand  and  capture  the  beam  coupling 
process.  In  this  paper  we  present  first  a 
comprehensive  review  of  neutralization  studies.  We 
then  present  a  series  of  simulations  using  a  2-D  PIC 
code1,2,3  as  well  as  the  implementation  of  a  3D 
PIC/DSMC  code.4,5  These  show  the  dependence  of 
the  beam  neutralization  on  beam  energy  and 
neutralization  current.  The  simulations  presented  in 
this  paper  serve  also  as  means  of  validation  of  the 
PIC-modules  of  our  PIC/DSMC  code  under 
development. 

HISTORY  OF  EP  NEUTRALIZATION 

Possibly  first  pointed  out  by  L.  Spitzer  in  1952 
[uncited  note  in  Seitz  et  al.  1961],  electric  propulsion 
plumes  needed  to  be  properly  mixed  with  electrons  or 
else  severe  space-charge  effects  would  result.  Before 
the  first  space  tests,  there  were  serious  doubts  as  to 
the  stability  of  any  neutralization  approach  to  the  ion 
beam  created  by  an  electrostatic  thruster.  The 
general  idea  was  for  neutralization  to  occur  shortly 
after  emission  to  prevent  beam  return.  However, 
there  was  lack  of  understanding  as  to  how  the 
electrons  would  stay  within  the  beam  if  they  were 
injected  or  if  the  neutralization  process  was  unstable 
to  small  perturbations.  Failing  to  properly  neutralize 
the  beam  would  cause  a  dramatic  reduction  in  thrust, 
as  a  significant  portion  of  the  beam  would  return  to 
the  spacecraft.  This  problem  was  first  addressed  by 
the  Ramo- Wooldridge  staff  in  their  review  of 
electrostatic  propulsion  in  I960.6  Their  one¬ 
dimensional  investigation  was  admittedly  unrealistic 
enough  to  provide  a  satisfactory  indication  as  to  the 
stability  and  practicality  of  neutralization. 

Looking  at  simplified  versions  of  the  problem, 
many  theorists  predicted  growing  instabilities  that 
could  turn  the  beam  back  to  the  spacecraft.  Seitz  et 
al.7  also  reviewed  the  process  in  1961,  noting  that,  if 
electrons  were  left  to  drift  into  the  exhaust  beam  from 
a  stationary  start  outside  of  it,  the  center  of  the  beam 
would  still  develop  a  virtual  cathode  and  suffer  from 
thrust  reduction.  On  the  other  hand,  they  found  that 
electrons  could  likely  oscillate  within  the  beam.  At 
the  same  time,  Baldwin8  pointed  out  the  necessity  of 
trapping  electrons  in  a  potential  well  created  by  the 


ion  beam  and  French9  found  that  oscillations  in 
electron  current  density  aided  in  neutralization. 
Pearlstein  et  al.  (1963)10  argued  that  there  would  exist 
a  possible  problem  if  the  neutralizer  could  only  emit 
in  a  space-charge  limited  regime  rather  than  an 
emission-limited  regime.  Mirels  (1961)11  argued  that 
the  location  of  emitter  and  thermal  motion  of 
electrons  would  not  significantly  impact  beam 
neutralization,  subject  to  more  general  design 
constraints.  Buneman  and  Kooyers12,  using  a  one¬ 
dimensional  PIC  code  in  1963  were  able  to  provide  a 
neutralized  beam  when  electrons  were  injected  at 
velocities  lower  than  the  directed  ion  velocity. 
Fluctuations  in  the  space  charge  field  provided 
mixing  of  the  beam.  Two  years  later  Wadhwa  et  al.13 
performed  a  two-dimensional  PIC  simulation 
showing  that  electrons  would  oscillate  within  the 
beam  to  allow  for  neutralization,  but  theorized  that 
oscillations  was  not  the  only  mechanism  at  work. 
One  method  suggested  was  that  fluctuations  in  the 
space-charge  field  allowed  for  entropy  increase  to 
mix  the  electrons,  but  these  fluctuations  were  not 
found  downstream  of  the  neutralizer. 

The  1964  Space  Electric  Rocket  Test  I  (SERT  I) 
found  that  it  was  quite  easy  to  neutralize  ion  beams  in 
space  from  a  simple  neutralizer  geometry.  In  a  series 
of  tests  it  was  shown  that  the  ion  thruster  developed 
thrust  at  a  level  indicating  complete  beam 
neutralization  even  with  the  neutralizing  filament 
biased  significantly  from  the  spacecraft  ground.14  In 
1970,  SERT  II  was  launched  with  the  objective  of 
demonstrating  long-term  ion  thruster  operation. 
While  it  quickly  developed  shorts  in  both  thrusters, 
these  were  resolved  by  1979  and  one  of  the  thrusters 
was  returned  to  operation.  In  reference  to  our 
discussion,  there  were  two  items  of  note:  The 
neutralizer  current  and  beam  current  were  nearly 
impossible  to  decouple  and  a  neutralizer  one  meter 
away  from  the  beam  was  able  to  neutralize  the  ion 
beam  indicating  that  neutralization  is  easier  than 
previously  thought. 

After  SERT  I,  proof  of  concept  was  achieved  and 
the  theoretical  discussion  of  beam  neutralization  was 
dropped  in  favor  of  engineering  new  thrusters. 
Studies  after  SERT  I  include  evaluations  of 
neutralizer  placement15'16,  optimization  of  the 
thrusters,  and  simulations  to  analyze  spacecraft- 
plume  interactions.17'18  Othmer  et  al.19'20  has 
investigated  the  neutralization  process  using  a 
relativistic  3D  PIC  simulation  to  model  neutralization 
of  particles  moving  at  around  0.1c,  some  three  to  four 
orders  of  magnitude  higher  than  typical  velocities. 
The  ratio  between  electron  and  ion  speeds  was  kept 
realistic,  however.  Their  findings  included  a  shock 
front  moving  through  the  beam  and  a  sheath  near  the 
neutralizer  injection  area  to  accelerate  the  electrons 
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into  the  beam.  Also,  a  strong  dependence  on  the  ratio 
of  electron  thermal  velocity  and  bulk  ion  speed  was 
observed  in  two  radically  different  modes  of 
neutralization.  The  relativistic  PIC  code,  as  well  as  a 
somewhat  unphysical  geometry  employed,  are 
fundamental  shortcomings  of  the  approach.  Tajmar 
and  Wang15  investigated  FEEP  neutralization.  They 
found  that  the  placement  of  the  neutralizer  had 
significant  impact  on  the  potential  distribution  of  the 
beam.  Further  numerical  studies  by  Tajmar  and 
Wang16  investigated  backflow  contamination  and 
charge  exchange  effects.  Neither  work  performed  a 
rigorous  study  of  the  electron  capture  process  and 
electron  motion  within  the  beam.  Similar  to  FEEPs 
are  Colloid  thrusters.  These  have  not  been  developed 
as  far  and  no  studies  of  colloid  neutralization 
simulations  exist,  to  the  knowledge  of  the  authors. 

Additional  flight  data  was  achieved  with  the 
Deep  Space  1  spacecraft,  which  recorded  over  16,000 
hours  of  use  of  its  NSTAR  Ion  Engine21.  Due  to  the 
fact  that  electrostatic  propulsion  was  a  proven 
concept,  less  interest  was  paid  to  neutralization  and 
most  numerical  studies  supporting  the  mission 
focused  on  contamination  and  sputtering  issues21  and 
in  fact  there  was  no  instrumentation  to  measure  the 
beam-neutralizer  coupling.22  Wang  et  al.18  used  a 
PIC  code  where  only  charge-exchange  particles  were 
modeled  as  particles  and  electrons  were  set  to  a 
Boltzmann  distribution  with  density  to  neutralize  the 
beam  appropriately.  In-flight  data  collection  focused 
on  charge-exchange  plasma  characterization  23 

Additionally,  ion  beam  release  studies  performed 
as  part  of  the  PORCUPINE  project24  suggested  that 
electrons  within  the  ion  beam  would  leave  after  a 
short  distance,  as  they  would  have  too  much  kinetic 
energy.  This  suggests  that  the  ambient  plasma  may 
play  a  greater  role  than  previously  assumed,  which 
may  change  operating  parameters  and  capabilities  of 
electric  propulsion  devices  in  different  environments. 

Despite  decades  of  research  and  the 
implementation  of  electric  propulsion  devices,  the 
detailed  process  by  which  an  ionized  beam  is 
neutralized  in  space  is  still  unknown.  Assorted 
methods  to  fit  data  with  theory  have  been  found,  but 
the  actual  process  has  yet  to  be  studied  in  sufficient 
detail  to  fully  understand  the  subject. 

UNSTRUCTURED  3D  PIC  CODE 
DESCRIPTION 

While  the  majority  of  our  simulations  were 
performed  using  the  code  XPDP21’2,3,  we  have  also 
begun  using  our  3D  PIC/DSMC  code  to  examine  the 
problem  in  a  more  realistic  fashion.  We  have 
developed  an  unstructured  grid  generator  that 
provides  three-dimensional  meshes  of  arbitrary 


Figure  1:  Tetrahedral  Cell 

geometry  and  allows  for  adaptation  of  the  mesh 
according  to  the  preliminary  solution  obtained  on  an 
initial  grid.4,5  The  generator  is  based  on  Watson’s25 
incremental  node  insertion  method,  using  properties 
of  Delaunay26  triangulation. 

The  general  procedures  for  loading  and  injection 
used  in  this  work  follow  Birdsall  et  al.27  and  Bird.28 

Integration  of  the  equations  of  motion  of  a 
charged  particle  are  performed  by  the  Boris  method29 
as  discussed  by  Birdsall  et  al.27  Particles  are  moved 
between  adjacent  cells  using  a  particle-tracing 
technique,  as  shown  in  Figure  1. 

To  formulate  a  finite  volume  method  for 
Poisson’s  equation 

Ni 

]T9in,  +  qenc 

=  =  - ,  (1) 

£»  £„ 

advantage  is  taken  of  the  Voronoi  dual  of  the 
Delaunay  triangulation  to  associate  an  irregular 
volume  to  each  node  on  the  grid.  The  Voronoi  cell 
corresponding  to  each  Delaunay  node  contains  the  set 
off  points  closer  to  that  node  than  any  other,  the 
facets  of  the  Voronoi  cell  are  orthogonal  to  the  lines 


Figure  2:  Delaunay- Voronoi  Cell 
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Constant  Simulation  Parameters 
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Simulation  Variables  and  Results 


67,1024 

7.6241 

111.8380 
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20.6984 

1.5144 


k=ihk 
A,  ... 


for  i=j, 


if  j  is  adjacent  to  i, 


7.67E+18 


7.67E+18 


1.21E+19 


1.4  1.21E+19 

1.5  3.83E+19  I  1.00 

1.6  |  3.83E+19  I  1.Q9  _ 

Table  1 :  Simulation  Parameters  for  Simulation  Set  1 

joining  the  tetrahedral  nodes  as  shown  in  Figure  2. 

This  method  reduces  Gauss’  law  for  a  node- 
centered  finite  volume  scheme  to  the  standard  2nd 
order  finite-difference  method  on  Cartesian  meshes. 

In  a  bounded  domain,  piece-wise  Dirichlet  and 
Neumann  boundary  conditions  specify  a  solution  of 
Poisson  s  equation.  Since  the  boundaries  of  the 
Delaunay  mesh  are  forced  to  coincide  with  the 
boundaries  of  the  computational  domain,  boundary 
condition  implementation  is  straightforward.  In  the 
case  of  a  Dirichlet  boundary  condition,  the  voltage  is 
placed  on  the  right  hand  side  of  the  matrix  and  the 
corresponding  row  is  zeroed,  with  a  one  placed  on 
the  diagonal.  Fluxes  in  Neumann  boundary 
conditions  are  added  to  the  flux  formulation  for  the 
Voronoi  cell  corresponding  to  the  boundaiy  node, 
with  the  value  of  the  inward  normal  electric  field 
multiplied  by  the  boundary  area  added  to  the  right 
hand  side  of  the  node  of  interest. 

In  matrix  form  with  boundary  conditions  as  in 
Figure  3,  Gauss’  law  is: 

f  1  0  0  0  HcM 


P 

T.21E-6 

1.21E-6 

4.85E-7 

4.85E-7 

4.85E-8 

4.85E-8 


Ri  j  =  0  otherwise. 

A 

~~  is  the  ratio  of  the  area  of 

the  Voronoi  face  Aq  between 
nodes  i  and  j  to  the  distance 
between  nodes  i  and  j  if  the 
nodes  L; The  boundary 


^2,1  ^2,2  ^2,3 

“^3,1  ^3,2  -^3,3 


^2, A r  4^2 

^3  ,N  •  ^3 


[Rn,1  Rn, 2  Rn,3  ”*  Rn,N  J 

Qi 

Q2  +  sqEN2AN2 
=  —  % 


N  is  the  number  of  nodes  in  the  mesh.  The 
coefficients  are  determined  by 


— I — 7-q5E-8 _ |  nodes  Ly.  The  boundary 

conditions  for  node  1  are 
Dirichlet  with  potential  <£>0 ,  and 
node  2  is  on  a  Neumann  boundary  with  inward  flux 

^N,2-^N,2  * 

SIMULATIONS  AND  RESULTS 
In  this  section,  we  present  the  results  of  a  series 
of  2D  and  3D  ion  beam  simulations.  As  explicit  3-D 
simulations  are  computationally  intensive,  the 
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Figure  3:  Boundary  Conditions  in  Delauny- 
Voronoi  Dual. 


Figure  4:  2-D  Simulation  Domain 
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majority  of  our  results  have  been  derived  using  the  2- 
D  code  XPDP2. 1,2,3 

2-D  Results 

We  performed  a  variety  of  simulation  runs 
examining  the  following  parameters:  current  density 
ji  and  je ,  injection  energy  EXi  and  EXe ,  and  thermal 

(random)  energy  Tt  and  Te .  The  cases  run  and  the 
selected  parameters  are  listed  in  Table  1.  The 
parameters  are  not  modeled  after  any  specific  system 
but  are  intended  to  be  representative  of  the  class  of 
problems  associated  with  electric  propulsion. 

The  simulation  domain  for  the  2-D  simulations  is 
shown  in  Figure  4.  A  square  domain  was  chosen  for 
simplicity  with  Lx  —  Ly  =  0.05  m  .  The  number  of 

cells  along  each  axis  was  ncx  =  ncy  =  128  .  A 

single  electrode  was  placed  along  the  left  wall  with 
length  Ldec  =  0.025  m .  A  background  density  of 

nbg  =  l.OxlO11  m'3  was  selected  to  compare  to 
LEO  operations.  The  particles  were  given  a 
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Figure  5:  Vx-x  phase  space  for  Case  1.1, 
t=1.8E-6s.  Top:  ions,  Bottom:  electrons. 


weighting  factor  of  1E4  to  allow  for  a  uniform 
background  of  simulated  particles. 

Results  from  the  2-D  runs  are  shown  in  figures  5 
through  11.  Figures  5  and  6  show  typical  Vx-x 

phase  space  plots  for  the  ions  of  a  neutralized  and  a 
non-neutralized  ion  beam  respectively.  Clearly 
visible  is  the  reflection  of  a  portion  of  the  beam  to  the 
emitter  in  the  non-neutralized  case.  The  slight  dip  in 
directed  velocity  as  visible  in  the  “neutralized”  case 
is  due  to  a  slight  potential  build-up  within  the  beam 
that  is  insufficient  to  retard  the  beam  to  reflection. 

First,  the  dependence  on  injected  electron  density 
was  examined  at  beam  energies  of  40,  100,  and  1000 
eV  and  a  plasma  density  of  1  x  1015  m‘3 .  This 
resulted  in  currents  in  the  milliamp  range,  i.e.  for  run 
1.1  the  ion  and  electron  current  was  matched  at 
qLdecLzj  —  3.07E-4  A  .  The  electron  current  was  set 
to  match  either  the  ion  beam  current  to  simulate 
current  neutralization,  or  so  that  the  resulting  density 
would  match  that  of  the  ion  beam  for  density- 
neutralization.  Throughout  this  series  of  simulations 
we  kept  the  electron  thermal  temperature  low,  but  not 


Figure  6:  Vx  -  x  phase  space  for  Case  1.2, 
t— 1.8E-6  s  .  Top:  ions.  Bottom:  electrons. 
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Figure  7:  Vx  -  x  phase  space  for  Case  1.3, 
t=1.8E-6  s .  Top:  ions,  bottom:  electrons. 
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Figure  8:  Vx  -  x  phase  space  for  Case  1.4, 
t=1.8E-6  s  .  Top:  ions,  bottom:  electrons. 


zero,  as  seen  in  Table  1.  Electron  directed  energy 
(  EXt )  was  set  to  zero  to  simply  “place”  the  electrons 
on  the  emission  surface. 

The  current-neutralized  beam  developed  a  virtual 
anode  and  reflected  at  40  and  100  eV  energies,  but  it 
was  neutralized  at  1000  eV,  showing  only  a  slight 
potential  increase  over  the  background.  Phase  space 
plots  are  shown  in  Figures  5  through  10.  The 
surprising  neutralization  at  1000  eV  beam  energy  led 
to  the  second  simulation  set. 

The  second  set  of  simulations  investigated  the 
relation  between  electron  temperature  and  beam 
energy.  A  requirement  for  low  electron  velocities  is 
similar  to  the  argument  put  forth  by  the  Ramo- 
Wooldridge  Staff.4  They  found,  using  a  one¬ 
dimensional  model,  that  the  beam  would  be  unstable 
unless  electron  speeds  were  less  than  ion  speeds.  A 
series  of  simulations  was  performed  to  find  where,  if 
anywhere,  the  electron  temperature  would  produce  a 
velocity  that  would  create  a  neutralized  beam,  in 
comparison  to  cases  1.1  and  1.3,  or  an  unneutralized 
beam,  in  comparison  to  case  1.5.  The  varied 
parameters  are  listed  in  Table  2,  all  other  parameters 
are  the  same  as  the  base  case. 

Unexpectedly,  the  lOOOeV  beam  would  only 
develop  a  virtual  cathode  when  the  electron 


temperature  was  brought  to  a  full  eV.  This 
corresponded  with  an  electron/ion  speed  ratio  of 
nearly  11.  The  lOOeV  case  could  not  become 
successfully  neutralized  even  with  an  electron/ion 
speed  ratio  of  0.5!  This  dramatically  illustrates  that 
electron  temperature  is  not  the  dominant  factor  in  the 
neutralization  of  ion  beams. 

In  all  simulations,  the  circulation  of  the  electrons 
along  the  beam  lengthwise  is  observed.  The 
velocities  of  the  circulating  electrons  are  significantly 
higher  if  they  are  unable  to  neutralize  the  space 
charge  of  the  ion  beam.  There  is  a  slower  “core” 
visible  that  corresponds  with  the  bulk  of  the  ion 
beam,  visible  in  the  unneutralized  simulations. 
Additionally,  one  can  see  a  sheath  formed  on  the 


Case 

Base 

case 

Te  (eV) 

V. 

(m/s) 

% 

Virtual 

Cathode? 

2.1 

1.5 

1.0 

419386 

10.94 

Yes 

2.2 

1.5 

0.1 

132621 

3.46 

No 

2.3 

1.3 

8.36E-4 

12123 

1.0 

Yes 

2.4 

1.3 

2.09E-4 

6061.5 

0.5 

Yes 

Table  2:  Parameters  and  Result  for  Simulation  Set  2. 
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Figure  9:  Vx-x  phase  space  for  Case  1.5. 
t=1.2E-6s.  Top:  ions,  Bottom:  electrons. 
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Figure  10:  Vx  -  x  phase  space  for  Case  1.6, 
r=1.2e-6  s  .  Top:  ions,  Bottom:  electrons 


emitter  surface  in  the  electron  phase  space  plots. 
Finally,  in  all  simulations  one  can  see  the  “surfing” 
effect  of  the  initial  ions  getting  accelerated.  This  was 
observed  by  Buneman  and  Kooyers10  as  well. 

As  can  be  seen  in  the  simulations  and  is  borne 
out  by  reality,  the  system  seeks  space  charge 
neutralization  at  the  cost  of  electron  momentum.  The 


Figure  11 :  Vx  -  x  phase  space  of  ions  for 
Case  2.3.  r=1.8E-6  s 


necessity  of  maintaining  a  sufficiently  small  excess 
of  electron  velocity  in  low  beam  energy  cases 
indicates  that  it  is  necessary  for  low-power 
neutralization  devices  to  operate  in  a  low-temperature 
regime,  permitting  the  electrons  sufficient  interaction 
time  to  slow  down  to  ion  velocities  before  exiting  the 
beam. 

In  all  cases  a  majority  of  the  electrons  do  get 
trapped  in  a  potential  well  created  by  the  ion  beam. 
Due  to  the  greater  energy  obtained  by  the  electrons  in 
unneutralized  cases,  the  “surfing”  acceleration  is  far 
more  noticeable  and  may  play  a  part  in  the  instability 
of  the  beam,  drawing  ions  farther  into  the  negative 
well  in  front  of  the  beam  and  negating  the  space 
charge  created  by  the  electrons  that  forced  them  to 
return  to  the  main  beam. 

The  role  of  the  perveance  P  in  beam 
neutralization  is  still  under  exploration.  While  it  is 
well  known  that  a  high  perveance  beam  will  develop 
a  virtual  cathode  much  sooner  than  one  with  a  lower 
perveance,  its  effect  on  the  ability  to  allow 
neutralization  has  not  been  fully  explored  yet.  Due  to 
the  dramatic  variance  in  susceptibility  to  electron 
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velocity,  perveance  likely  plays  a  large  role  in 
neutralization  stability. 

3-D  Results 

The  3-D  simulation  domain  is  shown  in  Figure 
12.  The  cylindrical  domain  of  radius  R  and  length 
L  consists  of  a  circular  emission  area  with  radius  r 
where  both  electrons  and  ions  are  injected.  While 
this  is  unphysical,  this  was  more  accurately  the  case 
performed  in  the  2-D  simulations.  Future  work  will 
include  a  separate  neutralizer. 

Due  to  the  longer  simulation  time  required  of  3- 
D  runs,  only  a  few  simulations  were  performed  for 
the  present  work.  A  domain  with  L  =  R  =  0.01  m 
was  generated  with  an  injection  surface  radius 
r  =  0.005  m .  The  background  was  held  at  a  density 
of  1E11  while  the  injected  density  was  1E15. 
Injected  velocity  was  set  to  10,000  m/s  for  the  first 

simulation  and  increased  to  100, 000  m/s  for  the 

second.  Injected  temperatures  were  held  at  0.1  eV 
for  both  ions  and  electrons. 


Figure  15:  3-D  Potential,  r=2.2E-7  s 


As  can  be  seen  from  Figures  1 3  and  1 6,  the  beam 
is  reflected  or  neutralized  under  the  same  conditions 
as  seen  in  the  2-D  case.  Again,  we  see  the  formation 
of  the  virtual  anode  reflecting  the  beam  when 
electron  velocities  are  higher  than  those  of  the  ions 
by  more  than  a  few  factors  of  two.  An  interesting 
result  is  in  Figure  15,  where  the  unneutralized  beam 
“arcs”  past  the  ions  before  returning  and  trapping 
very  efficiently  in  the  forward  region  of  the  ion 
beam.  This  indicates  that  the  electron  trapping 
mechanisms  are  of  significantly  different  strength  in 
3-D  as  compared  to  2-D. 

We  do  see  the  “surfing”  effect  observed  in  2-D 
here  as  well.  The  double  hump  of  the  accelerated 
head  corresponds  to  the  beginning  of  the  trapped 
electrons,  which  are  accelerating  the  head  of  the  ion 
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Figure  14:  3-D  Vx-x  phase  space  for  electrons. 
r=2.2E-7  s 
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beam.  In  the  Figure  14,  we  see  the  potential 
distribution  in  case  3.1.  The  virtual  cathode  is  plainly 
visible.  Its  location,  however,  is  strange  as  it 
corresponds  to  the  trapped  electron  concentration 
rather  than  the  location  of  the  slowest  ions. 

The  second  run  does  show  a  potential  well,  but 
the  ions  have  accelerated  out  of  it  without  the 
“surfing”  seen  in  the  first  run  indicating  a  more  stable 
solution.  Even  though  the  potential  shows  large 
increase  across  the  domain,  there  does  not  seem  to  be 
the  same  cathode  formation  as  we  have  seen  in  other 
cases. 


CONCLUSIONS 

The  history  of  electric  propulsion  neutralization 
development  was  presented  along  with  the 
outstanding  problems  therein.  Simulations  were 
performed  using  both  2-D  and  3-D  PIC  codes, 
demonstrating  computationally  the  neutralization  of 
an  ion  beam.  The  process  of  neutralization  was 
shown  to  be  dependent  on  the  electron  to  ion  velocity 
ratio.  This  compares  favorably  with  the  observed 
behavior  of  electrostatic  thrusters  and  the  theories 
presented  in  earlier  works. 
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Figure  16:  Vx  -  x  Phase  Space  for  Ions  in  3-D 
Case  2.  r=2.2E-7  s 


Figure  17:  Potential  contours  for  3-D  Case  2. 
r=2.2E-7  s 
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